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A dynamical model for quantum channel is introduced which allows one to pass continuously from 
the memoryless case to the case in which memory effects are present. The quantum and classical 
communication rates of the model are defined and explicit expression are provided in some limiting 
case. In this context we introduce noise attenuation strategies where part of the signals are sacrificed 
to modify the channel environment. The case of qubit channel with phase damping noise is analyzed 
in details. 
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I. INTRODUCTION 

In memoryless quantum channels successive signals 
(channel uses) are affected by independent, uniform 
sources of noise PISH SOI- O n ^ ne other hand, memory 
channels are characterized by the presence of correlated 
source of noise where each channel use is directly or indi- 
rectly affected by the previous ones. Preliminary results 
in the study of such systems has been obtained in Ref. (ij 
where it was pointed out that entangled codes can be 
useful in achieving optimal channel performances. Subse- 
quently some of these results have been generalized to the 
continuous variable case in Refs. 0,0, while a systematic 
analysis of the problem has been proposed in Refs. [9llTo|. 
In this paper we present a "dynamical" model for study- 
ing memory effects in quantum communication where the 
noise correlations are derived from the interactions be- 
tween the transmitted signals and the channel environ- 
ment. By varying the time intervals at which signals 
are produced by the sender of the message, the model 
simulates different communication scenarios. Memory- 
less configurations for instance are recovered as a limiting 
case in which the signals are transmitted at a frequency 
much lower than the inverse of the characteristic time 
of the channel environment relaxation. In this context 
we introduce also noise attenuation protocols where the 
sender alternates sequences of carrying-messages signals 
with sequences of signals which are employed to modify 
the environment response but which do not carry any 
messages to the receiver. Since timescales are fundamen- 
tal in our model, we characterize its efficiency by intro- 
ducing the transmission rates of the communication line. 
These are dimensional quantities (of dimension equal to 
an inverse time) which measure the maximal number of 
qubits or bits of information that can be transferred re- 
liably (i.e. with unit fidelity) through the channel per 
unit of transmitting time. Transmission rates are pecu- 
liar of our model as previous works |E S IE El were 
concerned in characterizing memory channels in terms 
of information capacities, i.e. the maximum number of 
qubits (or bits) that can be reliably transferred through 
the channel per channel uses. These figures of merit (i.e. 
rates and capacities) are in general distinct, but are pro- 
portional to each other when the sender of the message 



encodes her/his messages in regular sequence of signals 
(see Sec. lPTXL 

In Sec. [H] we introduce the channel model by focusing 
on the physical assumption which underline its definition. 
In Sec. IIIII we discuss the memory effects present in the 
system and we introduce the noise attenuation protocols. 
In Sec. IIVI and Sec. we define the transmission rates 
of the channel and we compute their values in some ex- 
tremal case. Finally in Sec. I VII an example of a dephasing 
qubit channel with memory is discussed. 



II. THE MODEL 

Consider a communication line where messages are en- 
coded into some internal degree of freedom (e.g. polar- 
ization, spin etc.) of a collection of identical physical ob- 
jects Ci, C2, • • • which propagate through the medium 
E that separates the sender (say Alice) from the receiver 
(Bob). The Cj are the information carriers of the sys- 
tem: they are locally produced by Alice and organized in 
a time-ordered sequence s = {n, T%, ■ ■ ■ } with Tj > be- 
ing the time interval between the instants tj+i and tj at 
which Cj+i and Cj enter E respectively. We will assume 
the effective transit time T tr it takes for the carriers for 
reaching Bob to be constant and shorter than the inter- 
vals Tj at which they are injected into the medium (fast 
propagation condition). The first condition guarantees 
that the time-ordering of s is preserved in the propaga- 
tion (i.e. Bob will receive the (J + l)-th carrier only after 
a time Tj from the arrival of the j-th carrier) . The sec- 
ond condition instead guarantees that E interacts only 
with one carrier at a time. Therefore, if R is the density 
matrix of the carriers at Alice location, after a time T tr 
Bob will receive the state 

R' = Ti E {W(R® Po ) W*} , (1) 

where po is the initial state of E, and where 

W = -~V j U j ■■■ V 2 U 2 V 1 U 1 , (2) 

is the unitary operator which describes the coupling be- 
tween the internal degree of freedom of the carriers and E. 
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FIG. 1: Schematic of the communication scenario. Alice en- 
codes her messages in the internal degree of freedom of the 
carriers Ci, C2, ■ ■ ■ , which propagates in a time-ordered se- 
quence toward Bob. The carriers interact one at a time with 
the local environment LE, while LE undergoes a dissipative 
evolution through its interaction with the reservoir R. 

In Eq. Q the terms 

IT* =Te>q>|-~ jf ' " dt [H CjE (t) + H E ]^ , (3) 

describe the interaction between Cj and E (here Hc^it) 
is the effective time dependent Hamiltonian that cou- 
ples Cj and E, while He is the free Hamiltonian of the 
medium). Working in a strong coupling regime we will 
neglect the contribution of He in Eq. J^J and we will 
assume the Uj to be uniform with respect to the label j. 
On the other hand, the terms Vj of Eq. (0) describe the 
free evolution of E in the time interval between the in- 
stant tj + T tr when Cj leaves the environment and the 
instant tj + i when Cj+i enters it, i.e. 

Vj = exp j - jH E {Tj - T tr ) j ^ exp j ~^H E Tj j . (4) 

In the following we identify two distinct components of 
the medium E: a finite dimensional Local Environment 
(LE) component which is directly coupled with the car- 
riers through the Uj, and a huge Reservoir (R) compo- 
nent which is coupled with LE but not with the carri- 
ers (see Fig. QJ. The free evolution 10} is supposed to 
induce a dissipative dynamics which transforms any ini- 
tial states of LE into a stationary configuration 00, with 
te being the characteristic time of the process. This 
is equivalent [TlT | to introducing a one-parameter family 
T = {£ t }t^o of Completely Positive, Trace preserving 
(CPT) which, given a the initial state of LE at some 
time to, represents its evolution at time to + r with the 
density matrix £ T (a). In this formalism £q coincides with 
identity map on TLle- On the other hand the stationary 
state Co of LE is defined by the property 

£ T (ao) = oo for all r ^ , (5) 




FIG. 2: Circuit representation of Eq. J7J. The local environ- 
ment LE interacts through the unitary couplings Uj (repre- 
sented by the small red circles in the figure) with one carriers 
at a time. Between two consecutive interactions with the car- 
riers instead LE undergoes the dissipative evolution described 
by the transformations £ T , (open circles). 



while the characteristic time te by the property 

S T>TE (e)=(T Tre, (6) 

for all bounded operator O of TLle- An example of T sat- 
isfying the above conditions will be presented in Sec. I VII 
Under the above approximations Eq. JQ) provides a 
bouncing ball description of the carrier-environment inter- 
actions where the carriers-balls move toward the LE-wall 
according to the time-ordered sequence s = {n, T2, • • • } 
chosen by the "pitcher" Alice and "hit" instantaneously 
the local environment LE one at a time (see Fig. . The 
resulting transformation is a time ordered product of in- 
teractions Uj and relaxation processes £ Tj (see Fig- EJ - 
Assuming LE to be initially in the stationary state (To 
this gives 

R' = Tr LB {■■■oS Tj oUjO... 

o£ T2 oU 2 o £ Tl 0U1 (i?(g)CT )} , (7) 

where the partial trace is performed on Hle, 
stands for the unitary mapping Uj(- ■ ■ )W on He, ^>H.le, 
and "o" indicates the composition of super-operators. It 
is important to note that in our model each sequence s = 
{ti, T2, • • • } is characterized by a distinct input-output 
relation Q. 



III. MEMORY EFFECTS 

Here we give an overview of the memory effects which 
are accounted for by the model introduced in Sec. ITT1 

Because of the time ordering of Eq. (0) the output 
state of a carrier might depend on the input state of the 
carriers which precedes it in s but it is always indepen- 
dent from the input state of the carriers which follows it 
in the sequence. As a matter of fact Eq. J7J) closely re- 
sembles the memory channels analyzed by Kretschmann 
and Werner in Ref. |10( . To make this more explicit we 
rewrite this equation in terms of of the discrete family of 
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CPT maps {$i n) }„ where 

$<?\R) = TTLE {Un°£r n ^°Un-l 

o-.-o£ Tl oWi(i?®<7 )}, (8) 

is the output state Q corresponding to the density ma- 
trix R of &j = iH.Cj associated with the first n carriers of 
the sequence s (here He is the Hilbert space associated 
with the internal degree of freedom of the j-th carrier). 
Therefore the model of Sec. HI originates prop er memory 
effects analogous to those of Refs. [ol IsLJa Ho| but avoids 
the feed-forward correlations of Ref. Q- For instance 
Markovian correlated noise can be recovered by properly 
choosing the transformations £ T (see Appendix lA"|l . 



A. Memoryless configuration 



Assume Alice is producing a sequence s with intervals 
Tj greater than or equal to the characteristic relaxation 
time te of the dissipation process T- see part a) of Fig. 
In this case, after each interaction, the local environment 
LE has enough time to relax into the stationary configu- 
ration Co before a new carrier begins interacting with it. 
Under this hypothesis Eqs. (JSJ) and (JSJ yield 



(9) 



where N is the CPT map which transforms the density 
matrices p of a single carrier into 



(10) 



Equation describes a memoryless configuration where 
the noise acts on the Cj independently. 



B. Generalized memoryless configuration 

A generalization of 10 is obtained when the carriers are 
organized in identical independent groups of m elements 
each. Here it is convenient to express the elements of 
s as T g j where g = 1,2, ••• is the group index, while 
£ € {1, • • • , m} labels the carriers within a given group. 
In this notation the time interval 



T 9,t 



(11) 



gives the "length" of the g-th. group while AT g — Tg :in 
is the interval which separates the last element of the g- 
th group from the first element of the (g + l)-th group. 
We do not assume any restrictions on the time inter- 
vals {T 9t i}i = i r .. , m _i which separates carriers belonging 
to the same group but we require carriers of distinct sub- 
groups to be separated by time intervals larger than te, 
i.e. AT g ^ te - see Fig. [3] part b). In this case from 
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FIG. 3: Some relevant configurations. Part a): memoryless 
configuration The carriers (represented by the green cir- 
cles) are separated by time intervals Tj which are greater than 
the dissipation time te of the local environment. Part b): 
generalized memoryless configuration 112|l . Here the carriers 
are divided in groups labeled by the index g. The groups are 
separated by time intervals AT g which are greater than the 
dissipation time te- Part c): perfect memory channel 1141 . 
Here the distance between two consecutive carriers is negligi- 
ble with respect to te inhibiting the relaxation of LE. Part d): 
example of a noise attenuation protocol. Alice sends uniform 
sequences of signals composed by n carriers (the B carriers 
of the protocol represented by yellow circles in the picture) 
which have been prepared in the same input state po and 
which are separated by time intervals r. These carriers do 
not convey any message to Bob and are employed only to 
"program" the environment response. The information is in- 
stead encoded into the (n + l)-th carrier (the A carriers of 
the protocol represented by the green circles). The sequence 
repeats after a time interval te to allow LE to return to the 
stationary configuration. 



Eq. (JSJ) follows that the transformation of the carriers of 
the first G groups can be expressed as 



$(») = 2>%xM 
where n = mG and 



(g) 



(12) 



M^{p)= (13) 
Tr LE {U g<m °^T S , m _i o ... o 8 Tg l oU gi i(p(g>a )} , 

is the CPT map associated with the m carriers C ffi i, 
of the g-th group. By comparison with 



C 



Eq. @, Eq. (|12|) describes a memoryless channel where 
the groups are the effective information carriers of the 
model. In particular if the sets {T gr i}e = t t ... jTra are uniform 
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with respect to the group label g, one has M. 



(ii) 



M 



0?') 



for all g and g' and the transformation (|12f) has once 
again the standard tensor structure Mf G . 



C. Perfect memory channel 

Consider the case where Tj <C te for all j. In this limit 
the local environment relaxation process is inhibited by 
the frequent interactions with the carriers. Consequently 
the £ Tj are replaced by the identity transformation on 
H LE and Eq. © yields 

$< B > (iZ) = Tr LE {U n o • • ■ o U 2 o U x (R <8> £T ) } • (14) 

This expression describes a perfect memory channel 0, 
Hfl | where the information transferred from the carriers to 
the finite dimensional local environment LE is not dissi- 
pated into the reservoir R of Fig. ^ These maps are 
asymptotically equivalent [lfj to noiseless channel where 
each carriers can transfer log 2 D qubits of quantum infor- 
mation reliably (here D is the dimension of the Hilbcrt 
space Tic 01 a single carrier). 



D. Noise attenuation protocols 

Here we present a communication strategy which ex- 
plicitly exploits the fact that in our model the environ- 
ment is effected by the signaling process. In this proto- 
cols only a subset A of the transmitted carriers is used to 
encode messages to Bob. The remaining carriers (subset 
B) are instead employed for perturbing LE in such a way 
that the Cj on which the messages are encoded have a 
better chance to reach Bob without being corrupted. In 
other words the B carriers are used by the sender as con- 
trol parameters to program the environment response to 
the A carriers. A simple implementation of a noise at- 
tenuation scheme is shown in Fig. [3] part d). Here the 
B carriers are composed by uniform strings of n states 
Po (represented by the yellow circles) separated by equal 
time intervals r. The information is instead encoded a 
single carrier (green circles) and the whole structure re- 
peats after a relaxation time te - this last assumption is 
not fundamental but allows us to treat the input-output 
relations of the A carriers as a memoryless channel of 
the form JJjJ. In this configuration the transformation of 
the A carriers which comes from solving Eq. Q can be 
computed as follows. First we determine the modified 
state o~ n of LE which arises from the interactions with 
the B carriers. This is accomplished by solving the set of 
coupled equations analogous to those of Ref. [13] , 



a'j = Tr c {U(p (So*)}, 
o-j+i = £ r (cr'.) , 



(15) 



where the trace is performed over the carrier degree 
of freedom, 11 is the usual carrier-LE coupling super- 
operator and j = 0, 1, • • • , n — 1. The density matrix 



a n which results from i|15|) is then used to determine the 
output state of the A carriers according to the equation 



W{p) = Tv LE {U(p(g>a n )} . 



(16) 



The transformation (|16|) is in general different from 
Eq. (|10|) and depends explicitly on the parameters n, 
r and po that are controlled by Alice. The basic idea 
of a noise attenuation scheme is to appropriately select 
such parameters in order to get a transformed mapping 
A/" which is less noisy than the original mapping J\f. An 
example of this effect will be presented in Sec. IVII 



IV. TRANSMISSION RATE OF A SEQUENCE 

Timescales play a fundamental role in the model pre- 
sented in Sec. [H] Therefore a proper way to characterize 
it, is by introducing its quantum and classical transmis- 
sion rates. In simple terms these quantities measure, re- 
spectively, the maximum number of qubits and bits per 
second that Alice can encode into the carriers sequence 
s without compromising the readability of the transmit- 
ted messages. The formal definition of the rate of the 
sequence s is constructed as follows. 

First of all we introduce the discrete value function 
n s (T) which, given the sequence s, counts the number 
of carriers which fit 13] in the time interval [0, T[. Fur- 
thermore, for any e > and T > we define q s (e,T) to 
be the dimension -in qubits units- of the largest Hilbcrt 
sub-space of TL(T) = 8^1^ which allows for a fi- 
delity of the transmitted state greater than 1 — e. This 



q s (e, T) = max \ log 2 d : 3H co de dimTt co de = d, 3 A, V 

d V. 

V|tf) eH code F($,Do$f)oi)>i- e }, (17) 

where TL co de are Hilbert sub-spaces of Tt(T), A and D 
are encoding and decoding CPT maps on Ti(T) applied, 
respecitively, by Alice and Bob to the carriers, and 

F{^,Vo&p oA) = (*|Do$P oA(\y}{y\)\V) ,(18) 

is the fidelity between the input state |\&) 6 TL co de and 
the decoded output state Vo $i T ^ oA{\^)(^\) (for easy of 
notation $i T ^ indicates the map 

$ K(T)) of E ^ p that 

acts on the n s (T) carriers of s which lie on [0, T[). The 
quantum transmission rate r q (s) of s is thus given by the 
ratio q s (e, T)/T in the the limits e — > 0, T — > 00 , i.e. [3 



j(s) — lim limsup 



qs(e,T) 
T 



(19) 



Analogously we define the classical transmission rate 
r c (s) of s by substituting the function q s (e, T) with 
the largest number of classical distinguishable messages 
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c s (e,T) that can be transmitted to Bob with fidelity 
greater than 1 — e, i.e. 



lim limsup 



r c (s) 

where as in Eq. (|17fl one has 



c 3 (e,T) 
T 



(20) 



c s (e, T) = max \ log 2 d : 3H co de dimH co de = d, 3 A, V 

~ik £{!,■■■ ,d} F(* fe ,Do$p oA) > 1 - e }(21) 
with |^i), |^2), • • • j \^d) being an orthonormal basis of 



first-neighbors distance among the carriers. These are 
the regular sequences of the model: for them Eqs. i|22|) 
and l|25|) coincide and the transmission rate is propor- 
tional to the quantum capacity of the channel, i.e. 



r q (s) = Q s /t s 



(27) 



The same analysis can be repeated also for the classi- 
cal rate r c (s) of Eq. (|20|l . In particular, in this case, 
Eqs. J22J, an d 123 still apply by replacing Q s with 
the classical capacity C s of the maps {$i }„ defined by 



C' 



lim lim sup 

e >0 n — >oo 



c s (e,n) 



(28) 



A. Upper and lower bounds 

A simple upper bound for the quantum rate r q (s) of s 
can be derived from Eq. i|19[l as follows, |15| 



r q (s) = lim limsup 



n s (T) q s {e,T) 
T n s (T) 



lim lim sup - 4 lim sup 

e^O t^cc n s (T) 1 T '->oo 

Q.K > 



n s (T) 
T' 



(22) 



where is the minimum average first-neighbors distance 
among the carriers of s defined by 



■tii v n s (T') n s (t) 
1/t s = limsup — — — = hm sup 

T'->oo T T'^oc t ^ T , t 



On the other hand 



n v v Qs{e,T) q s (e,n) 

Q s = lim lim sup ; — — = lim limsup (24) 

e^O T^oo n s (T) c-*0 n^oo n 

defines the quantum capacity 0, El E, associated 
with the maps } n of Eq. JSJ (in this expression 

q s (e,n) is given by l(T7|) with n s (T) replaced by n). 

A lower bound for r q (s) is instead obtained as fol- 
lows fil 



r q (s) — lim limsup 



n s {T) g s (e,T) 
T n s (T) 
q s {e,T) 



lim limsup 



lim inf — - 

T'^oo T' 



(25) 



where r" ^ is the maximum first-neighbors average 
distance among the carriers of s defined by 

1/t'J = lim inf = lim inf ^ . (26) 

If the sequences s is such that hin^^oo n s (T)/T = l/r s 
exists, one has = r" = r s with r s being the average 



B. Some solvable configurations 

The maximizations implicit in Eqs. (|24|) and (j28(l are 
in general difficult to solve. However, following the anal- 
ysis of Refs. 0, 0] one can bound the capacities Q s 
and C s by means of the coherent information 0] and 

of the Holevo information of $1 , respectively. In 
particular we have 



lim sup max 



J& N \R) 
N 



(29) 



where the maximization is performed over all density ma- 
trices R of N carriers and 



(23) J(^ N) ,R) = S($i N \R)) - ®I A )(*r)), (30) 



is the coherent information 0] of (i?). In the above 
expression = — Tr[i?log 2 i?] is the von Neumann 

entropy, *S?r is a generic purification of i? constructed 
by adding an ancillary Hilbert space TIa, and Xa is the 
identical map on TLa- Analogously one has 



C s ^ limsup max 

N^oc V 



x(^ N \v) 

N 



(31) 



where the maximization is performed over all ensemble 
= {Pk] Rk}k of N carriers and where 

7>) = S($MC£ Pk R k )) (32) 

k 

-J2PkS($i N \R k )) , 

k 

is the Holevo information associated with $ S N \ 
Kretschmann and Werner have identified a class of maps 
{$s }„ (the forgetful channels f° r which the right- 
hand side term of and (|3*T)l indeed provide the exact 
value for Q s and C s . Here we will focus only on the limit- 
ing cases discussed in Sec. lIIII for which an expression for 
Q s and CL can be derived without the elegant arguments 
of Ref. [lfl. 
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a) The simplest configuration is when the sequence s 
is such that Tj <C te for all j. When this hap- 
pens the maps {<&i"^}„ describe a perfect memory 
channel (|14fl which allows optimal transfer, ensur- 
ing Q s — C s — log 2 D. Therefore, according to ljT7|) 
using regular sequences s with Tj <C te, Alice and 
Bob can achieve transmission rates equal to 



r q (s) = r c (s) 



\og 2 D 



(33) 



b) For memoryless configurations @, Q s and C s co- 
incide, respectively, with the quantum Q(Af) and 
classical C (AT) capacity of the memoryless map AT 
of Eq. I|10|) . On one hand one has Q , 



Q(AT) — lim max 

N^oo R 



J(Af® N ,R) 
N 



(34) 



where, as in Eq. 1)29(1 the maximization is per- 
formed over all density matrices R of N carriers and 
where J(Af® N , R) is the coherent information 130|) 
f On the other hand one has H, 



C(Af) = lim max 



N 



(35) 



where the maximization is performed over all en- 
semble V = {pk',Rk}k of A carriers and where 
x(Af® N ,V) is the Holevo information (|32[) associ- 
ated with Af® N ■ Therefore for regular sequences s 
with Tj ^ te we get 



r,(a) = QCAO/7 



r c (a) = C(AA)/7 



(36) 



The generalized memoryless configurations l|12JI 
can be treated in the same way by replacing the 
quantities t' s ,t" of Eqs. and (J2SJ) with the 

corresponding average first-neighboring group dis- 
tances and the map M with the m carriers memory- 
less map A4 S of Eq. i|13|) . In particular, for a gen- 
eralized memoryless sequences s having constant 
group lengths T g = T s and constant group separa- 
tions AT g = AT S for all g one easily verifies the 
following identities 



r,(a) = Q{M S )/(T S + AT S 
r c (s) = C(M S )/(T S + AT S 



(37) 
(38) 



d) Finally consider the noise attenuation protocols of 
Sec. IIII Dl For the sake of simplicity we will focus 
on the specific example of Fig. where the results 
for memoryless configuration applies. In this case 
the rate is given by 



r q( S ) 

r c {s) 



Q{N)/{nT- 
C(N)/{nT - 



te) , 



(39) 



with AT being the map (|16f) and with tit + te being 
the time intervals which separates two consecutive 
A-carriers. 



V. TRANSMISSION RATE FOR MULTIPLE 
CHOICE OF THE SEQUENCE 



In this section we analyze the optimal quantum and 
classical communication rates i? 9jC achievable in our 
model when Alice is not restricted to a single given se- 
quence s but instead she has some freedom in selecting 
the sequence she will use for the signaling. 

For the sake of simplicity we will assume the set S 
of the allowed sequences to be fully characterized by a 
single parameter r m i„ which bounds the minimum value 
for the intervals Tj of a sequence s of the set. That is 
S = S(r m i n ) will be the set of all sequences s which sat- 
isfy Tj T m i n for all j. The need of constraining the min- 
imum value of the Tj is fundamental if we want our model 



to have a non trivial structure (see for instance Sec. lIVBl 
and Eq. I|46l) below). From a more practical point of 
view the introduction of T m i n follows from the physical 
and technological difficulties in producing sequence of or- 
dered signals that might arise in realistic communication 
scenarios (for instance, too close packed carriers tend to 
overlap during their propagation, compromising the time 
ordering of the sequence) . 

A natural candidate for R q , c is the maximum of the 
rates r qtC (s) computed over the sequence s of S, i.e. 

R^liTmm) = max r q , c (s) ■ (40) 

A detailed analysis of R q ]l is presented in Appendix [Bl 
where it is shown how Eq. (|40|l simplifies in the case 
in which S contains only regular sequences for which 
Eq. (|27|l applies. We will see in a moment that for 
Tmin <C t e and 

Train 

^ t e the function i%i(r m , n ) pro- 



vides indeed the correct values of the achievable rates. 
For generic r m , n however we claim that the function 
Rq]l{T m i n ) does not necessarily tell the whole story about 
R q ,c- On the contrary we propose to compute R q , c as fol- 
lows 



R„ 



= 



lim lim sup max 



Rc{T m in) = hm lim sup max 



ses 



qs(e,T) 
T 

Cs(e,T) 
T 



(41) 
(42) 



with q s (e,T) and c s (e,T) given in Eqs. (TJJ and JUJ. 
Equations l|41|) and (|42|l define proper rates of the com- 
munication line of Sec. [I] in the sense that, given <5 > 
and e arbitrarily small there is allowed sequence s G S 
which, in the limit of infinite T permit Alice to transfer 
to Bob at least (R q — S)T qubits with fidelity > 1 — e. 

Since Eq. igDJl is obtained from Eqs. igTJ an d 63 by 
inverting the order of the maximization over s with the 
limits in e and T it follows immediately that R q ]l(T m i n ) 



is a lower bound for R qiC {T m in) of 5, i.e. 

R q ,c(Tmin) ^ Rq.c^Trnin) • 



(43) 



An interesting problem is to understand whether or not 
the inequality in Eq. (|43|l can always be replaced with an 
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identity. Alternatively one may ask under which condi- 
tions on the model parameters (i.e. Uj, T) the transmis- 
sion rate of S can be computed as the maximum of the 
rates achievable within a specific choice of s. In the next 
section we provide a partial answer to these questions by 
showing that for Tmin t e and r m i n ^ te the functions 
Rq,c{r m in) and -Rg!c(r mm ) coincide. 



A. Bounds and asymptotic behavior 



Even without solving the maximizations of 14UI) . (|41|> 



R 



Q,C 



and (|4*2|) one expects the resulting expressions R q ]l 
will depend upon the interplay between the relaxation 
time te of LE and the characteristic time T m ; n of S. 

A trivial but useful upper bound for R qc follows by 
observing that the maximum number n s (T) of carriers 
that can fit in [0, T[ cannot be greater than T/T m i n and 
that q s (e,T),c s (e,T) cannot be greater than the log 2 of 
the dimension of Tt(T), i.e. 



q s (e,T),c s (e,T) ^n s (T)log 2 D, 



(44) 



with D being the dimension of the Hilbert space of a 
single carrier. Replacing the above relations in Eqs. I|19|) 
and H2Ufl gives 



Rq ,c(^"mm) ^ 



log 2 L> 



(45) 



for all 



i . From Sec. HVBl it follows that this bound 
is achievable at least if S is such that T m i n « T£. In 
this case in fact the sequence Sq with Tj — r min for all j 
allows for carriers that reliably transfer log 2 D qubits of 
information each. Therefore from (|40|l and (|43() we get 



H (t ■ 1 = R^(t ■ ) 

- LL q,c\ 1 min J — Jt g,cv 1 min ) 



\og 2 D 



,«T E 



(46) 



which shows that the rates diverge for r m i„ — * 0. An ex- 
plicit expression can also be determined for T m j„ greater 
than te- In fact, according to Sec. IIII Al in this case 
all the allowed sequences s yields the same memoryless 
mapping J\f® n ( T \ Thus the maximization with respect 
to s becomes a simple optimization with respect to the 
average time intervals r s and one gets, 



Rqi^~min) — Rq {^~rriin) 

R ( T ■ 1 = ffWfr ■ "1 

J - L c\'min) ±L c \'minj 



Tmin ^T E 



Q(A0/r mm ,(47) 
C(A0/r min ,(48) 



with Q(Af) and C(N) the capacities of Eqs. (53} and (|o3jl . 
respectively. 

For intermediate value of r m .;„ a lower bound for R q ]l , 
and thus for i? 9iC , can be obtained for instance by fo- 
cusing on the generalized memoryless configuration (see 
Eq. I|B1|) ) or by considering the noise attenuation strate- 
gies. In this last case it is simpler to consider only the 
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FIG. 4: Plot of the ratio F of Eq. I|56|l as a function of the 
dimensionaless parameter t/te, for different values of the n 
and for different values of the environment-carriers coupling 
constant A. In the strong coupling regime A ~ 0, the atten- 
uation noise protocol provides a significative improvement of 
the transmission rate. For instance for A = 0.01, r reaches 
the maximum value of ~ 1.3 for n — 1 and r ~ te/2. 



configurations described in Fig. [3] and maximizing the 
rates l|39[) with respect to the free parameters r ^ T m i n 
and n ^ 1, e.g. 



q yT mm ) ^ SUp 



i?«( 



Q(F) 



T^T min nr + t e 

n>l 



> sup 



C(N) 



r^T min nr + te 



(49) 



VI. AN EXAMPLE WITH QUBITS 

In this section we analyze an example of dynamical 
model for memory channels where both the information 
carriers Cj and the local environment LE are qubits. In 
this context we will make a comparison between the noise 
attenuation protocol of Sec. IIII Dl and the memoryless 
configuration. 

We will assume the carrier-LE interaction Uj of Eq. Q 
to be to a control-unitary such that when the carrier is in 
|0)c7 nothing happens to LE, while when Cj is in |l)c 
the environment undergoes to the transformation 



©(A) 



\/T — A 

s/T^X -Va 



(50) 



with A G [0, 1] being a parameter which measures the 
"intensity" of the coupling (with low coupling corre- 
sponding to A ~ 1 and high coupling corresponding to 
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A ~ 0). Moreover we will assume the relaxation process 
T = {£ T } T acting on LE to be described by amplitude 
damping maps 1] which takes the state \1)le to \0}le 
with probability 1 — t)(t) where rj(r) € [0, 1] is a non 
increasing function of r with characteristic time te, i.e. 



£A\0)le(0\) = \0)le(0\ 

£ t (\1)le(1\) = v{t)\1)le(1\ + 0-v(t)) \0) LE ( 
£A\0)le(1\) = VW)\0)le(1\ ■ 



(51) 



In this example the stationary state <jq of LE is hence 
|0)lb. The parameterization of the memory effect is 
given by rj(r), with r\ = corresponding to the memory- 
less case (fast environment relaxation) and rj = 1 corre- 
sponding to perfect memory case (no environment relax- 
ation). In order to have a well defined threshold between 
memoryless and memory configuration, in the following 
we will assume 



1 — t/te for t < te 
for r ^ t e . 



(52) 



Under the above conditions, it is possible to show that 
both the map N of the memoryless case and the map 
J7 l|16fl of the noise attenuation protocol correspond to a 
phase damping channel V g where the coherence terms of 
the input qubit p are degraded by a positive factor g ^ 1 , 



7VI«)c<«|) 

r g (\o) c (i\) 



l«)c(«| 
9 |0)c(l| 



for k = 0, 1 



(53) 



In particular Eq. I|10|) gives Af — V go with go = On 
the other hand, Eq. Ijl6|) gives Af — Vg where ~g is a com- 
plicated expression (jClOfl of A and of the parameters po, 
n and r (see Appendix IO for details). By appropriately 
selecting the values of the above quantities one can make 
the make Af less noisy than Af by having ~g > go . To see if 
this corresponds to an increase in the transmission rates 
r q,c{s) we can use the results of Sec. IIVBI In the case 
of the phase damping channels V g the capacities Q(V g ) 
and C(V g ) of Eqs. ffity and l|3"5|) can be explicitly com- 
puted. For instance since here the noise does not affect 
the populations associated with the computational basis, 
the classical capacity of the phase damping channel Ij53(l 
is optimal for all values of g, i.e. C(V g ) = 1. Hence from 
Eqs. © and we get 



; (so) = 1/te ^ l/(nr + te) = r c 



(54) 



where so is the memoryless sequence with uniform inter- 
val Tj = te and r c is the classical rate of the noise at- 
tenuation protocol of Fig. [3] Equation H54|) shows that, 
in the specific example considered here, the noise atten- 
uation protocol does not improve the classical rate of 
the communication line with respect to the memoryless 
case. On the other hand the quantum capacity of a phase 
damping channel l|53[) is equal to [2(j 



Q(V g ) = l-H 2 (l/2 + g/2) 



(55) 



where H 2 (x) = — xlog 2 x — (1 — x) log 2 (l — x) is the bi- 
nary entropy function. In this case, higher values of g 
corresponds to higher Q(V g ) and the rate r q of the noise 
attenuation protocol can be higher than the rate r q (so) 
of the memoryless case. To see this we studied the ratio 



r q = T E l~ff 2 (l/2 + g/2) 
r q (s ) ut + t e I- H 2 (l/2 + go/2) ' 



(56) 



as a function of the variable t/t e and for for different 
values of n and A. [Here g has been optimized with re- 
spect to the no-carrying signal po] ■ The results have been 
plotted in Fig. 0] which shows that in the strong coupling 
limit A ~ one can have an appreciable increase of T for 
t ~ te/2 and with n of the order of 5. 



VII. CONCLUSION 

We have introduced a communication model where 
memory effects arise from the interaction between the in- 
formation carriers with the channel environment. Differ- 
ent memory effects can be simulated by varying the time 
intervals at which the carriers are produced by the sender 
of the message. The information rates of the model have 
been defined and computed in some extremal cases. 



APPENDIX A 



In this appendix we show how a Markovian correlated 
noise 0, 0, [HJ can be derived from the mapping (JHJ by 
properly choosing the transformation £ T . . 

Consider the case in which for sufficiently big r the 
map £ T describes a decoherent process of LE where, given 
an orthonormal basis airihE^ one has 



£t(\1)le{Z'\) =5e,e> |^(r))z B <^WI 



(Al) 



with the vectors {| - 0^( r )) le}i being not necessarily or- 
thogonal, and 8i£i being the Kronecker delta. The con- 
dition (JHJ can then be satisfied by identifying cfq with 
one element of the selected basis (say \Iq)le), and im- 
posing | ipg (t ^ te)) le = \£o) le for all I. In this case the 
mapping (JHJ can be expressed in terms of the operators 



A,, = 



LE 



LE 



{hW^io) 



LE 



(A2) 



(£ j+1 \U j+1 \^(Tj)) LE , (A3) 



which act, respectively, on the Hilbert space rtc-i ano - 
fic j+1 for j = 1, • • • , n — 1. They allow us to define the 
probability distribution 

p£> = Tt Ci {4A (i } (A4) 

and the conditional probabilities 

4t'k = Trc J+1 {4 +1 ,,A + ^} ■ ( A5 ) 
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which satisfies the normalization conditions 
£/ p¥ +1 l = 1 and J2f P^ +1 l < !■ Using 
these quantities Eq. ijHJl can be finally expressed in 
compact Markovian form, 

*M ln<tn _ x ■ ■ ■ M e2/l M tx R M\ y ■ ■ ■ A// n _ Mn 

with Mt t = I \J P^ti an< ^ 

Me J+1 ^=A 4j+1 ,t j / y /^^ j . 

APPENDIX B 

In this section we analyze R q ]l showing that, if the set 
S contains only regular sequences, then the maximization 
of Eq. H4l)f) can be solved by focusing on the generalized 
memoryless configurations. 

Consider the subset Sq of the sequence s £ S which 
correspond to the uniform generalized memoryless con- 
figurations of Sec. IIII 51 characterized by constant group 
distance AT S = max{T m i„, t e }. Since Sq is a proper 
subset of S we have 



R { q\T m in) > max r q (s) 



ses 



Q(M S ) 



ses (> T s + max{r mm , t e } 



(Bl) 



where we used Eq. (|37f) to express r q (s). Now, given 
s e S from Eqs. lt2*2l and one gets 



r q {s) < (l/T^)limsupmax(j($( JV) , J R)/iv} 

J([<t> { s N) }® k ,R') 



^ (1/t' s ) lim sup < sup max • 

A<"->oo k>l Rl 



(l/r' s ) hmsup < lim max 

N->oz |fc->oo R' 

(l/r s ) hmsup ^ 



J([$ 



kN 



kN 



(B2) 



where in the second and in the third line the maximiza- 
tion is performed over the density matrix R' of k x N 



carriers, [$i Af - ) ]® fc are k copies of the map &s"\ and 
Q($s ) is the memoryless quantum capacity (|34[l of the 
map $s . The second inequality is trivial: it follows 
from the fact that maxjj J(&i N \ R)/N can be obtained 
from max ff J([$£ JV) ]® fe , R')/(kN) for k = 1. The iden- 
tity on the third line instead is a consequence of the fact 
that maxff J([<f>i N) ]® k ,R')/(kN) achieves its maximum 
for k — > oo. We can further simplify the above expression 
by introducing the time interval T S (N — 1) = X^jLi 1 T j 



(AO 



associated with the first N — 1 carriers of the sequence s 
and noticing that 



,. T S (N-1) „ 
hmsup = t s 



(B3) 



with r" defined as in Eq. (|26|l . Using this result, from 
Eq. ijEl we get 

T s (iV - 1) + max{r mi „, t e } 



r q (s) ^ lim sup 



x lim sup — — 

at^oo T S (N - 1) + max{ T min ,T E } 

^ t^. N T S (N - 1) + max{T mm ,r E } 



sup 

-7"' .- ,-, 



t£ se5o T s + max{T min ,T E } 



(B4) 



The ratio t' s '/t' s is always greater than or equal to one. 
However, if the set S includes only sequences which are 
regular, than for all s we have t' s = r'J. In this case the 
bounds of Eqs. IjBljl and l)B4|l coincides yielding 



Q(M t 



ses ( , T s + max{r mra , t e } 



(B5) 



The same derivation applies also for the classical rate 
i?c^. In this case one can show that if S contains only 
regular sequence then, 



i? (1) (Y ■ ) - 



max 



C(M S 



gSo T s + max{r mi „, t e } 



(B6) 



a. Asymptotic limit 

It is interesting to note that the above expressions give 
the correct asymptotic values of Sec. IV Al For instance 
for T min ^ t e we have A4 S — J\f® m where m is the num- 
ber of carriers contained in each group of the sequence 
and TV is the memoryless map @. Given s 6 5o this 
yields 



mQ(JV) 

T s + max{r min , t e } T s + 



Q(M S ) 



where we used the additivity property Q(AT® m ) = 
rnQ(Af) of memoryless channels and the fact that group 
length Ulljl is always greater or equal to (m — l)r m i n . 
Equation 1471) finally follows by noticing that the rate 
Q(.N~)/T m i n is achieved by the sequence of iSo with r g ^ — 
Train for all <? and I. 

The limit (|46|l instead follows by noticing that the 
rate log 2 £>/r m j n can be obtained from the set Sq by 



using T g 



for all I 



of large group, i.e. 



= 1, • • • , m — 1 in the limit 
oo. In this case in fact 
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M s is a tensor product of perfect memory channels and 
T s = (m - l)r mm , so that 



Q(M S 



mlog 2 D 



T s +max{T mm ,T £ } 



(m - l)r min + r B 
log 2 £> 



(B8) 



APPENDIX C 



To characterize the modified map of J7 we first solve 
the system <|15f) by using the following parameterization 
for the density matrices element of <jj in the canonical 
basis {\0}le, \1)le], 



°J X 3 + l Vj 

iy 3 zj 



(CI) 



with Zj G [0, 1] and Xj, yj real for all j — 0, 1, • • • ,n. 
The resulting recursive equation can be simplified by in- 
troducing the column vectors 



w = (1 - p){^\\ - A), Tj^^MT^Xjf 
and the 2x2 Hermitian matrix 



A=(l-p) 



r,(^ - 1 + 2A) -2ry 3 / 4 yA(T 



A) 

1-2A) 



where 77 stands for rj (r) and p is the population associ- 
ated with the \0}c component of the no-carrying message 
state po. In this notation Eq. (fT5|l gives the following un- 
coupled equations 



Vj+i = Vrj ( 2 P ~ !) 



A • Vj + w 



(C2) 
(C3) 



which can be solved analytically. In particular, imposing 
the initial condition ctq = |0)lb(0| (i.e. xq = ya = zq — 
0) the first one gives yj — for all j. The solution of (|C3|) 
instead can be obtained in terms of the eigenvalues A± of 
A and their corresponding eigenvectors (a±, P±) T '. Ex- 
plicitly the eigenvalues of A are 



Ah 



Vv 



[(1 + Vv)p 



-(l-p)(l-07)(l-2A)±A] 



(C4) 



with 

A = {4V»7(l-2p) (C5) 
+ [(1 + v^)p + (1 -P)(l - Vv)(l ~ 2A)] 2 } 1/2 • 

The corresponding eigenvectors (a±,/3±) have instead 
the following components 



a± = r t 1 l\l-p)^\{l-\)/N ± , (C6) 

ftfc = [(V^-l)p-(l-p)(l-2A)(l + V^)=FA]/JV±, 
with the normalization coefficient 

Ar ± = {16(l-p) 2 V^A(l-A) (C7) 
+ - 1)P - (! - P)(l - 2A)(1 + v/A) t A] 2 }V2 . 
In particular, for |A±| < 1 one has [2^J| 



3-1 



fc=0 



1 - A> 
1-A 



(C8) 



and thus 



Zj = ?7 3/4 (l-p) r ? 1 / 4 (l-A)u^') + v /A(T — A)i W 



xj = 7? 1/2 (l-p) ry 1/4 (l-A)^''') + v / A(l-A)z; (j) 



where vP> = « a) = C+Vif + f-^-, and 

t (3) = fk') a+/9+ + e 0) a_/3_ with 



1-A± 



(C9) 



Setting j — n and replacing the above expressions into 
(|C1() we obtain the modified state of LE a n after n suc- 
cessive interactions with pg. Using the definition Ijltifl one 
verifies that AT is a phase damping channel (|53(l charac- 
terized by a damping factor 



VA - 2 (VA z n — Vl — A x n ) 
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